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Abstract. Hypertree decompositions of hypergraphs are a generaliza- 
tion of tree decompositions of graphs. The corresponding hypertree- 
width is a measure for the cyclicity and therefore tractability of the 
encoded computation problem. Many NP-hard decision and computa- 
tion problems are known to be tractable on instances whose structure 
corresponds to hypergraphs of bounded hypertree- width. Intuitively, the 
smaller the hypertree- width, the faster the computation problem can be 
solved. In this paper, we present the new backtracking-based algorithm 
det-fc-decomp for computing hypertree decompositions of small width. 
Our benchmark evaluations have shown that det-fc-decomp significantly 
outperforms opt-fc-decomp, the only exact hypertree decomposition algo- 
rithm so far. Even compared to the best heuristic algorithm, we obtained 
competitive results as long as the hypergraphs are not too large. 



1 Introduction 

Since many important problems in computer science are intractable in general, 
it is a competitive task to identify tractable subclasses of such problems which 
can be solved efficiently. One approach to do this is to restrict the structure of 
a problem represented as graph or hypergraph. For example, the structure of 
instances of the constraint satisfaction problem CSP and the equivalent Boolean 
conjunctive query problem BCQ can be naturally encoded by hypergraphs. A 
hypergraph is a generalization of a graph where each edge connects a set of 
vertices. Gottlob et al. [3] have shown that in analogy to tree-width of graphs, 
the hypertree-width of hypergraphs is an appropriate measure for the cyclicity 
and therefore the tractability of the corresponding computation problems. The 
hypertree-width of a hypergraph is defined as the minimum width over all de- 
compositions of the hypergraph in a way called hypertree decomposition. Roughly 
speaking, the smaller the width of a hypertree decomposition, the faster the cor- 
responding problem can be solved. However, deciding whether there exists a 
hypertree decomposition of width at most k is NP-complete in general. 
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Gottlob et al. [3] have shown that for fixed k, the problem of deciding whether 
there exists a hypcrtrcc decomposition of width at most k is in LogCFL, i.e., 
it can be solved in polynomial time and is highly parallclizable. Moreover, they 
presented the alternating algorithm fc-decomp, which constructs a hypertree de- 
composition of minimal width less than or equal to k (if such a decomposition 
exists). Another algorithm for computing hypertree decompositions of minimal 
width less than or equal to k is opt-fc-decomp [2, 7, 9]. Up to now, opt-fc-decomp 
has been the only implementable exact polynomial-time algorithm for construct- 
ing fc-bounded hypertree decompositions. However, opt-fc-decomp has an im- 
portant disadvantage: although polynomial, it needs a huge amount of memory 
and time even for small hypergraphs. This basic problem remains even after im- 
provements like redundancy elimination as investigated by Harvey and Ghose [5] . 
Thus, recent research focuses on heuristic approaches for constructing hypertree 
decompositions of small but not necessarily minimal width. One of the practically 
most successful algorithms of this kind is due to McMahan [8], who combined 
well-known tree decomposition and set cover heuristics in order to construct hy- 
pertree decompositions. This approach, however, has the disadvantage that the 
computation cannot be focused on hypertrees of width smaller than some upper 
bound k, i.e., the algorithm returns a hypertree decomposition whose width may 
be much larger than the minimal one although there is often time left to improve 
this result. 

In this paper, we consider a combination of exact and heuristic approaches 
in the sense that we restrict the search space by a fixed upper bound k and 
apply heuristics to accelerate the search for a hypertree decomposition of width 
at most k (but not necessarily the minimal one). Such an approach has the fol- 
lowing advantages: the search space and the width of the resulting hypertree 
decomposition can be restricted by k. In particular, this means that the algo- 
rithm is able to find a hypertree decomposition of minimal width by setting k 
small enough. On the other hand, if k is not minimal, the algorithm may con- 
struct a hypertree decomposition of non-minimal width but within a reasonable 
amount of time since the solution space becomes larger. In principle, such an 
approach allows us to control the tradeoff between the resulting hypertree width 
and the required computation time. Thus, it combines the advantages of exact 
and heuristic algorithms while it minimizes their disadvantages. 

We implemented our algorithm det-fc-decomp based on these insights. Our 
experimental evaluations have shown that it performs much better than opt-fc- 
decomp. In particular, it needs only a small amount of memory and is much 
faster than opt-fc-decomp. Moreover, it is often the case that det-fc-decomp 
finds hypertree decompositions of width smaller than those obtained by McMa- 
han's heuristic algorithm. Finally, there is a further advantage of our algorithm: 
due to its top-down nature of decomposing a hypergraph into sub-hypergraphs 
and decomposing the sub-hypergraphs into sub-sub-hypergraphs etc. it can in 
principle be implemented for parallel execution in order to increase its perfor- 
mance even further. 
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Fig. 1. Example of a hypergraph and its hypertree decomposition of width 2 



This paper is organized as follows: In Section 2, we give the basic defini- 
tions used in this paper. Afterwards, in Section 3, we describe the alternat- 
ing algorithm fc-decomp. Then, in Section 4, we introduce our new algorithm 
det-fc-decomp and prove its correctness and polynomial runtime. In Section 5, 
we present our evaluation results. Finally, we conclude in Section 6. 

2 Preliminaries 

A hypergraph H is a tuple (V, E), where V is a set of vertices (variables) and 
E C 2 V \{0} is a set of hyperedges (atoms); w.l.o.g., we assume that V\U E = 0. 
We define vertices(H) = V and edges(H) = E. A hypertree for a hypergraph H 
is a triple (T,x,A), where T = (V,E) is a tree and X ■ V -> 2 vertlces(H ^ and 
A : V — > 2 edaes ( H ^ are labeling functions. We define vertices (T) = V and refer 
to the vertices of T as "nodes" to avoid confusion with the vertices of H. If 
V = (V',E') is a subtree of T, we define x(T') = \J peV , xip)- We denote the 
root of T by root{T), and for every p 6 vertices (T) we denote the subtree of T 
rooted at p by T p . 

A hypertree decomposition of a hypergraph iJ is a hypertree (T, x, A) for H 
satisfying the following four conditions [3]: 

1. Ve G edges(H) 3p G vertices(T) : e C x(p) 5 

2. Vw G vertices(H) : the set {p G vertices(T) | w G x(p)} induces a (connected) 
subtree of T, 

3. Vp G vertices (T) : x(f) ^ (JA(p), and 

4. Vp G vertices(T) : [J \(p) n x(T p ) C x(p)- 

The OT<i</i o/a hypertree decomposition (T, x, A) is given by max pet , el , ttces (y) |A(p)|, 
and the hypertree-width of a hypergraph is the minimum width over all its hy- 
pertree decompositions. Fig. 1 shows an example of a hypergraph (on the left) 
and its hypertree decomposition of width 2 (on the right). 

Let H = (V, E) be a hypergraph. A path between two vertices x, y G V 
is a sequence of vertices x — v\, t>2, ■ • ■ , Vk-it v k — y such that for all pairs 
Vi,Vi+i there exists e G E with v i} v i+ i G e. Now, let W C V. A set V C V is 
[VF] -connected if for all x, y G V there exists a path between x and y not going 
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through vertices in W. A [W]- component is a maximal [VF]-connected non-empty 
set of vertices V C V \ W. We call the set W in this context a separator. Note 
that both, components and separators, can also be represented by the hyperedges 
covering the corresponding vertices, which we will do in this paper. 

3 The Algorithm fc-decomp 

In this section, we describe the alternating algorithm fc-decomp introduced by 
Gottlob et al. [3]. This algorithm checks non-deterministically whether there ex- 
ists a hypertree decomposition of a hypergraph of width at most fc. By using 
this algorithm, Gottlob et al. proved that the problem of deciding whether a hy- 
pergraph has fc-bounded hypertree- width is in LogCFL, a very low polynomial- 
time complexity class which is contained in NC2 and which thus consists of 
highly parallclizablc problems. The algorithm works in a top-down manner by 
guessing in each step the A- labels and checking two conditions. The x-labels 
and sub-hypcrgraphs can then be computed deterministically. A game-theoretic 
intuition for such an approach can also be found in [4]. 

The algorithm fc-decomp shown in Algorithms 1 and 2 is a notational vari- 
ant of the algorithm presented in [3] . A proof of equivalence of the two variants 
is given in the appendix. The main procedure fc-decomp in Algorithm 1 ex- 
pects a hypergraph HGraph consisting of a set of vertices vertices (HGraph) 
and a set of hyperedges edges (HGraph) as parameter. It calls the recursive 
procedure fc-decomposable in Algorithm 2, which returns a hypertree decom- 
position of width at most k if it exists and NULL otherwise. The parameters of 
fc-decomposable are a set OldSep of hyperedges which were chosen as separator 
in the previous run and a set Edges of hyperedges which represents a component 
w.r.t. DldSep, i.e., which represents the current sub-hypergraph that has to be 
decomposed. In line 1, a set Separator of hyperedges of size at most k is guessed. 
In lines 2 to 4 it is then checked whether Separator is indeed a separator, i.e., 
whether it decomposes the current sub-hypergraph into several (at least one) 
smaller components and whether the conditions of a hypertree decomposition 
are satisfied. In line 6, the procedure separate computes these components and 
stores them in Components. Afterwards, in lines 7 to 15, the hypertree decom- 
positions of the components are recursively computed and stored in Subtrees. 
Line 16 computes the x-labels of the current hypertree node in such a way that 
the conditions of a hypertree decomposition are satisfied. Finally, the procedure 
getHTNode in line 17 constructs the current hypertree node consisting of the 
A-labels represented by Separator, the x-labels represented by Chi, and the 
subtrees stored in Subtrees. 

For example, when applying fc-decomp to the hypergraph in Fig. 1, fc-de- 
composable is called with parameters Edges = {A, B, C, D, E, F, G, H} and 
DldSep = 0. In the first run, we guess Separator = {A, B} which results in 
Components = {{C, D, E}, {F, G, H}}. Consequently, there are two recursive 
calls: the first one uses parameters Edges = {C,D,E} and OldSep = {A,B}, 
and the second one uses Edges = {F,G,H} and OldSep = {A, B}. In the fol- 
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lowing recursion level, we guess Separator = {C, D} and Separator = {F, H} 
respectively. It is easy to see that no further recursive calls are necessary since 
both sub-hypergraphs are then completely decomposed. The resulting hypertree 
decomposition is shown in Fig. 1. 



Algorithm 1 k-decomp(H Graph) 

1 HTree := k-decomposable(edges(HGraph),$); 

2 return HTree; 



Algorithm 2 k-decomposable{Edges 1 OldSep) 

1 guess Separator C edges(HGraph) such that | Separator\ < k; 

2 check that the following two conditions hold: 

3 |J Edges n (J OldSep C \J Separator; 

4 Separator n Edges 7^ 0; 

5 if one of these checks fails then return NULL; 
e Components := separate(Edges, Separator); 

7 Subtrees := 0; 

a for each Comp G Components do 

9 HTree := k-decomposable(Comp, Separator); 

10 if HTree = NULL then 

11 return NULL; 

12 else 

13 Subtrees := Subtrees U {HTree}; 

14 endif 

15 endfor 

16 CTw := (|J Edges n (J OldSep) U |J(5'eparator l~l Edges); 

17 HTree := getHTNode(Separator , Chi, Subtrees); 
is return HTree; 



As it was shown in [3], /c-decomp can be implemented on a logspace ATM 
having polynomially bounded tree-size. Hence, deciding whether a hypergraph 
has fc-bounded hypertree- width is in LogCFL and thus in P. However, since 
fc-decomp is a non-deterministic algorithm, this is only of theoretical interest 
and cannot be implemented. The only deterministic algorithm with polynomial 
runtime for computing fc-bounded hypertree decompositions published so far is 
opt-fc-decomp [2, 7, 9]. This algorithm, however, has a basic drawback: due to its 
bottom-up approach, it needs a huge amount of memory even for small examples. 
So it is practically only applicable to hypergraphs of very small size and simple 
structure. Our aim in the following section is therefore to transform fc-decomp 
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Fig. 2. Dependency graph between procedures called by fc-decomp 



into a deterministic algorithm with polynomial runtime which performs better 
than opt-fc-decomp. We call our resulting algorithm det-fc-decomp. 



4 The Algorithm det-fc-decomp 

In this section, we present our new algorithm det-fc-decomp which we obtain 
from fc-decomp by replacing the guess and check-part in lines 1 to 4 of Al- 
gorithm 2 by a backtracking-based search procedure. In order to preserve the 
polynomial runtime, however, we have to store already visited separators; other- 
wise, the same subtree could be constructed multiple times which would lead to 
an exponential runtime. Algorithm 3 represents the main procedure which calls 
the recursive procedures decompCov, decompAdd, and decompSub described in 
Algorithms 4 to 6. We divided the recursive part into these three sub-procedures 
for better readability. Their dependencies are shown in Fig. 2, where an arrow 
from A to B means that A calls B. Moreover, we use three auxiliary procedures: 

— separate is the same as for fc-decomp described in Section 3 and computes 
the components of the current sub-hypergraph w.r.t. the chosen separator. 

— cover selects hyperedges for the separator which are satisfying the first con- 
dition of the check-step in Algorithm 2; we will describe it later in more 
detail. 

— expand is used to complete a "pruned" hypertree decomposition: if the same 
hypertree node is constructed the second time during our search procedure, 
we know already if the corresponding components can be decomposed. So we 
have to prune the hypertree at this node in order to guarantee the polynomial 
runtime as mentioned above. The procedure expand expands such nodes to 
subtrees after the search process, which can be done in polynomial time. 

We are now going to explain the transformation from fc-decomp to det-fc- 
decomp in more detail. The main procedures in Algorithm 1 and Algorithm 3 are 
almost the same. The only differences are the instantiation of the sets FailSeps 
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Algorithm 3 det-k-decomp(HGraph) 



1 FailS eps := 0; 

2 SuccSeps := 0; 

3 HTree := decompCov(edges(HGraph),%); 

4 if HTree NULL then 

5 HTree := expand (HTree); 
e endif 

7 return HTree; 



Algorithm 4 decompCov (Edges, Conn) 

1 if \Edges\ < k then 

2 //Tree := getHTNode(Edges, (J Edges, 0); 

3 return HTree; 

4 endif 

5 BoundEdges := {e 6 edges(HGraph) \ e n Conn ^ 0}; 

6 for each CovSep £ cover (Conn, BoundEdges) do 

7 HTree := decompAdd(Edges , Conn, CovSep); 
s if Tree / iVt/XL then 

9 return HTree; 

10 endif 

11 endfor 

12 return NULL; 



Algorithm 5 decompAdd(Edges, Conn, CovSep) 

1 InCovSep :— CovSep l~l Edges; 

2 if InCovSep / or fc - CowSepI > then 



3 if InCovSep = then AddSize := 1 else AddSize := endif; 

4 for each AddSep C Edges s.t. | J 4dciS'ep| = AddSize do 

5 Separator := CovSep U AddSep; 

6 Components := separate(Edges, Separator); 

7 if VComp £ Components . {Separator, Comp) ^ FailSeps then 

8 Subtrees := decompSub (Components, Separator); 

9 if Subtrees ^ then 

10 CTii := Conn U U(^ re C 0? ''S' e P U AddSep); 

11 HTree := getHTNode(Separator , Chi, Subtrees); 

12 return HTree; 

13 endif 

14 endif 

15 endfor 



16 endif 

17 return NULL; 
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Algorithm 6 decompSub(Components, Separator) 



1 Subtrees := 0; 

2 for each Comp G Components do 



3 ChildConn := \J Comp n |J Separator; 

4 if (Separator , Comp) € SuccSeps then 

s //Tree := getHTNode(Comp, ChildConn,®); 

6 else 

7 HTree := decompCov(Comp, ChildConn); 
g if FTree = NULL then 

9 FailSeps := FailSeps U {(Separator , Comp)}; 

10 return 0; 

11 else 

12 SuccSeps :— SuccSeps U {(Separator , Comp}}; 

13 endif 

14 endif 

15 Subtrees := Subtrees U {//Tree}; 



16 endfor 

17 return Subtrees; 



and SuccSeps, which are used to store already visited separators, and the call 
of expand. If a component could be successfully decomposed, the corresponding 
separator together with the component is inserted into SuccSeps; otherwise it 
is inserted into FailSeps. If a new separator is constructed during the search, it 
is checked whether it occurs in one of these two sets before its components are 
tried to be decomposed. This avoids that the same components are decomposed 
multiple times. Intuitively, since there are at most 0(n ) separators each of size 
at most k and each with at most n — \edges(HGraph)\ components, the search 
procedure runs in polynomial time for fixed k. 

Now, let us consider the recursive part in Algorithms 4 to 6. First note 
that the second parameter of decompCov is a set of vertices instead of a set 
of hyperedges in fc-decomposable. This is just a simplification since OldSep in 
Algorithm 2 occurs only in the expression (J Edges n (J OldSep. So we compute 
this set of "connecting vertices" Conn before the recursive call and replace OldSep 
by Conn. 

The first few lines in Algorithm 4 are a simple optimization: if the current 
component contains at most k hyperedges, it can be trivially decomposed into a 
single hypertree-node. Afterwards, in line 5, our determinization of the guess and 
c/iecfc-part in Algorithm 2 starts. In particular, we compute the set BoundEdges 
of hyperedges which are necessary to satisfy the condition in line 3 of Algo- 
rithm 2, which is Conn C (J Separator. This condition can only be satisfied 
by the separator through hyperedges in BoundEdges, i.e., through hyperedges 
containing some vertices in Conn. Thus, we can reduce the search space for an 
appropriate separator by first selecting a subset CovSep of at most k hyperedges 
in BoundEdges such that Conn C (J CovSep, where CovSep C Separator. This 
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selection is done by the procedure cover in line 6. In particular, cover succes- 
sively returns all possible selections and we have to loop through all of them. 
The further decomposition steps are then performed by decompAdd, which is 
called in line 7. If one such decomposition is successful, the resulting hypertree 
is returned in line 9; otherwise, NULL is returned in line 12 after all selections for 
CovSep have been tried. Note that the order of choosing CovSep is crucial for the 
duration of the search process in the case of success. Thus, we use heuristics in 
cover to obtain an appropriate selection order. We will explain these heuristics 
later in more detail. 

Let us now consider the procedure decompAdd defined by Algorithm 5. Its 
third parameter is the set CovSep which contains at most k hyperedges. These 
hyperedges were selected in such a way that the condition in line 3 of Algo- 
rithm 2 is satisfied. Since every separator guessed in fc-decomposable must sat- 
isfy this condition and we loop through all possibilities of CovSep, we know that 
our search space is large enough to find all separators which are accepted by 
fc-decomposable. This, however, is not necessary. To improve the runtime of 
our algorithm, we can restrict the search space based on the following insights. 

Definition 1. [3] A hypertree decomposition of a hypergraph is in normal form 
if for each r G vertices (T) and for each child s of r all the following conditions 
hold: 

1. there is (exactly) one [x{r)]~ component C r such that x(T s ) — C r U (x(s) fl 



2. x( s ) nC r ^ 0, where C r is the [x{ r )]~ component satisfying Condition 1, and 



Definition 2. Let (T, x, A) be a hypertree decomposition of a hypergraph and 
s G vertices(T). We say s is minimally labeled if s = root(T) and \X(s)\ = 1 or 
r is the parent of s (with component C r ) and it holds that for each e G A(s): 



Definition 2 means that removing any edge from the A-labels violates the 
connectivity or the monotonicity condition. 

Definition 3. A hypertree decomposition of a hypergraph is in strong normal 
form if it is in normal form and each vertex is minimally labeled. 

Lemma 1. For each k-width hypertree decomposition of a hypergraph H there 
exists a k-width hypertree decomposition of H in strong normal form. 

Proof. First note that we can assume w.l.o.g. that our hypertree decomposi- 
tion (T, x, A) is in normal form [3]. We will show the lemma by an inductive 
argument on the size of the x-components. To this aim, let s G vertices(T) be 
not minimally labeled and assume that all other vertices on the path from s 
to root(T) are minimally labeled. Now, we distinguish between two cases as il- 
lustrated in Fig. 3: (a) If s = root(T), we create a new vertex r' and add it to T 



X(r)) 



3- lU(s)n X (r) C X ( S ). 



1. U edges(C r ) f~l U A(r) % \}{\{s) \ {e}); or 

2. (A(s) \ {e}) n edges(C r ) = 0. 



( Connectivity) 
(Monotonicity) 
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(b) 







Fig. 3. Transformation steps to strong normal form 



as parent of s. We set A(r') = {e} for any hyperedge e £ A(s). (b) Otherwise, if 
s ^ root(T), let r be the parent of s. We create a new vertex r' and add it to T 
between s and r, that is, r' is the parent of s and r is the parent of r'. We set 
A(r') = A(s) and remove hyperedges from A(r') until any further removal would 
violate one of the conditions in Definition 2. Finally, we set x( r ') = x(s)fllj A(r') 
and remove all v e x(r) \ x(r') from x( s ') i0r eaC11 s ' G vertices(T s ). It is easy 
to see that the resulting tree is still a hypertree decomposition, however, it is 
not necessarily in normal form. Since our modifications affect only the subtree 
rooted at r', we know that this subtree can be transformed into normal form 
according to [3] without changing the labels of r' . Since C r i C C r by the mono- 
tonicity condition, our inductive argument applies. □ 

This means that it suffices to consider the smallest separators satisfying the 
conditions in Algorithm 2. In particular, for decompAdd in Algorithm 5, this 
implies that it suffices to add at most one hyperedge to CovSep such that the 
condition in line 4 of Algorithm 2 is satisfied. To this aim, decompAdd computes 
in line 1 the intersection between CovSep and Edges. In line 2 it is then checked 
if the condition in line 4 of Algorithm 2 can be satisfied at all. In particular, if 
InCovSep is empty then we must add an edge to satisfy the condition. However, 
if CovSep contains already k hyperedges, this is not possible and we have to reject 
CovSep. Otherwise, in line 3 we set AddSize to 1 if wc must add an edge and 
to if CovSep already satisfies the condition. Then in line 4 we loop through all 
hyperedges in Edges if AddSize = 1, that is, if we must add an edge; otherwise, 
AddSep is the empty set and the loop body is executed only once. Finally, our 
separator is computed in line 5 as the union of CovSep and AddSep, where AddSep 
contains either a single hyperedge or is empty. In the latter case, CovSep satisfies 
both conditions in Algorithm 2 (recall that it satisfies the first condition by 
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construction). Based on Lemma 1, we know that our computation of Separator 
is a correct determinization of the guess and c/iecfc-part in Algorithm 2. 

The call of separate in line 6 of Algorithms 2 and 5 is now completely 
the same in both algorithms. Since each sub-component returned by separate 
must be decomposable, we check in line 7 of Algorithm 5 if one of these sub- 
components is already known to be undecomposable. If so, we abort and choose 
an alternative separator. Otherwise, we try to decompose the sub-components, 
which is done in the procedure decompSub in line 8. The resulting hypertree 
decompositions are returned in Subtrees if the decomposition was successful 
for all sub-components; otherwise, Subtrees = 0. Note that there must be at 
least one sub-component, i.e., Subtrees cannot be empty in the case of success. 
Otherwise, a trivial decomposition of the current component would have been 
computed in lines 1 to 4 of Algorithm 4. The construction of a new hypertree 
node is then performed in lines 10 and 11 completely analogous to lines 16 and 17 
of Algorithm 2. 

Finally, let us consider the procedure decompSub defined in Algorithm 6, 
which contains the functionality of lines 8 to 15 of Algorithm 2. In line 2 we loop 
through all components that have to be decomposed and in line 3 we compute 
the set of vertices which will be the second parameter Conn of the recursive call 
of decompCov. In order to guarantee the polynomial runtime, we have then to 
check in line 4 if the current sub-component is already known to be decompos- 
able. If so, we prune the hypertree in line 5; otherwise, we apply the recursive 
call in line 7 analogous to line 9 in Algorithm 2. Depending on the result, i.e., 
whether the sub-component could be decomposed or not, we store it either in 
FailSeps or in SuccSeps. Moreover, the constructed hypertree decompositions 
are stored in Subtrees in line 15 analogous to line 13 in Algorithm 2. 

It is thus easy to see that our algorithm det-fc-decomp is a deterministic 
variant of fc-decomp defined in Algorithms 1 and 2. This insight of the correctness 
of det-fc-decomp is now formulated in the following: 

Lemma 2. For any given hypergraph H such that hw(H) < k, det-fc-decomp 

accepts H. Moreover, for any c< k, each c-width hypertree- decomposition of H 
in strong normal form is equal to some witness tree for H. 

Proof. Let HD = (T, x, A) be a c-width hypertree decomposition of a hyper- 
graph H in strong normal form, where c < fc. We show that there exists 
an accepting computation tree for det-fc-decomp on input H which coincides 
with HD. To this aim, note that by Lemma 9 in [3] there exists an accepting com- 
putation tree for fc-decomp on input H which coincides with HD. Thus, it suffices 
to show that every accepting computation tree for fc-decomp with corresponding 
hypertree decomposition in strong normal form is also an accepting computa- 
tion tree for det-fc-decomp. This, however, is trivial since det-fc-decomp tries all 
computation trees of fc-decomp whose corresponding hypertree decomposition is 
in strong normal form. □ 
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Lemma 3. If det-fc-decomp accepts a hypergraph H, then hw(H) < fc. More- 
over, each witness tree for H is ac-width hypertree- decomposition of H in normal 
form, where c < k. 

Proof. We show that every accepting computation tree for det-fc-decomp on 
hypergraph H coincides with some c-width hypertree decomposition of H in 
normal form, where c < fc. To this aim, note that by Lemma 13 in [3] every 
accepting computation tree for fc-decomp on input H coincides with some c- 
width hypertree decomposition of H in normal form, where c < fc. Thus, it 
suffices to show that every accepting computation tree for det-fc-decomp is 
also an accepting computation tree for fc-decomp. This, however, is trivial since 
det-fc-decomp tries a subset of all computation trees of fc-decomp. □ 

By combining Lemma 2 and Lemma 3 we get: 

Theorem 1. det-fc-decomp accepts a hypergraph H if and only if hw{H) < fc. 

Another important property beside the correctness of det-fc-decomp is its 
runtime. We gave already some hints why it runs in polynomial time, but we 
will now consider this in more detail. For simplicity, we assume in analogy to [7] 
that the standard set operations can be performed in constant time. Moreover, 
for better comparability, we ignore the auxiliary procedure expand since it is not 
necessary for the decision problem. 

Theorem 2. det-fc-decomp runs in time 0{n k+1 min(n, ck) k to 2 ), where n is 
the number of hyperedges in the hypergraph and c is the maximum number of 
incident hyperedges over all hyperedges. 

Proof. First note that there are at most 



/V / \ A, 



^ i\(n — i)\ 



separators each with at most to sub-components. Thus, since we restrict the 
number of recursive calls by checking SuccSeps and FailSeps to the number of 
possible separators and sub-components, the number of recursive calls is bounded 
by 0{&rri). What remains to show is the runtime of a single recursive call. To 
this aim, note that the number of loops in line 6 of Algorithm 4 is bounded by 

k k 

/min(n, cfc)\ min(n, ck)\ 

~ ^ I i ) r-i i!(min(n,cfc) - i)\ 

Moreover, the number of loops in line 4 of Algorithm 5 is bounded by n and the 
number of loops in line 2 of Algorithm 6 is bounded by m. Thus, the runtime of a 
single recursive call is bounded by 0(<Pnm). The overall complexity is therefore 
0(\P <Pnm 2 ), which in turn is bounded by 0(n k+1 min(n, ck) k to 2 ). □ 



12 




(a) opt-fc-decomp (b) det-fc-decomp 



Fig. 4. CPU and memory usage when decomposing s298 of the ISCAS89 benchmarks 
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Fig. 5. Covering Conn = {a, b, e} by BoundEdges = {^4, B, D, F, G, H} 



Since we can assume that cfc is in many cases much smaller than n, we have 
shown that det-fc-decomp has a better runtime complexity than opt-fc-decomp, 
which runs in time 0(n 2k m 2 ) [7]. Moreover, note that similar to opt-fc-decomp, 
the runtime may be significantly smaller in practice. Another important fea- 
ture of det-fc-decomp is its low memory usage compared to opt-fc-decomp. For 
example, the CPU and memory usage of both algorithms when applied to exam- 
ple s298 of the ISCAS89 benchmark suite can be seen in the screenshot of the 
Microsoft Windows Task Manager in Fig. 4. Example s298 describes a simple 
circuit hypergraph consisting of 133 hyperedges and 139 vertices. In Fig. 4 (a), 
we see that the memory usage suddenly increases to its maximum when ap- 
plying opt-fc-decomp; therefore, the system resources are busy with memory 
swapping such that there are only minimal resources available for solving our 
decomposition problem. After an hour opt-fc-decomp was still trying to solve 
the problem without success. In contrast in Fig. 4 (b), there is an unrecognizably 
small usage of memory when applying det-fc-decomp; therefore, the full CPU 
power can be used for our decomposition problem, which was successfully solved 
within 90 seconds. We present more experimental results in Section 5, which 
demonstrate that det-fc-decomp significantly outperforms opt-fc-decomp. 
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decompCov({A, B, C, D, E, F, G, H}, 0) 
decompAdd({A, B, C, D, E, F, G, H}, 0, 0) 
decompSub({{B, C, D, E, F, G, H}}, {A}) 
decompCov({B, C, D, E, F, G, H}, {a, b, c}) 
decompAdd({B , C, D, E, F, G, H}, {a, b, c}, {A}) 
decompSub({{C, D, E}, {F, G, H}}, {A, B}) 
decompCov({C, D, E}, {a, c, d, /}) 
decompAdd({C, D, E}, {a, c, d, }}, {C, D}) 
decompSub({{E}} , {C, D}) 
decompCov({E}, {g, i}) 
decompCov({F, G, H}, {a, b, e}) 
decompAdd({F, G, H} , {a, b, e}, {A, F}) 
decompSub({{G, H}}, {A, F}) 
decompCov({G, H}, {a, e, h}) 



A = {A} 
■ = {a,b,c 



A = {A, B} 
-- {a, b, c, d, e, /} 



{o 



{C.-D} 
C, d, f, g 



= {E} 



A = {A,F} 
; = { a, b ; e, h } 



A = {G, H} 
C = {a.,e,h,j} 



Fig. 6. A procedure call sequence and the corresponding hypertree decomposition 



Now, let us come back to our heuristics in the auxiliary procedure cover. As 
we mentioned before, the order of choosing CovSep in line 6 of Algorithm 4 has a 
crucial influence on the performance of det-/c-decomp. We obtained our results 
with a very simple heuristic approach exemplified in Fig. 5 by a decomposition 
step when decomposing the hypergraph in Fig. 1. Let us assume that Conn is 
given by {a, b, e} and BoundEdges is given by {A, B, D, F, G, H}. The task of 
procedure cover is now to compute subsets of BoundEdges that cover Conn such 
that at least all minimal covers arc considered. Recall that the minimal covers 
are sufficient because of Lemma 1. We do this in the following way: at first we 
assign to each hyperedge in BoundEdges a weight which is the number of vertices 
in Conn it contains. This can be seen on the left of Fig. 5 where the weight of A 
is 2 since A n Conn = {a,b}, the weight of B is 1 since B n Conn = {e}, etc. 
Then we order the hyperedges according to their weight as it can be seen in 
the last line on the left of Fig. 5. Finally, we choose the separators based on 
this ordering in the following way: we consider the hyperedges from left to right 
and select a hyperedge which covers a vertex in Conn that is not covered by 
other hyperedges selected so far. On the right in Fig. 5, the selection of the five 
separators {A,F}, {A, B}, {A, G}, {F, D}, and {F,H} is shown. For example, 
the first one is selected by choosing A which covers a and b and F which covers e; 
the second one is selected by choosing B instead of F to cover e; etc. All five 
separators are computed in this way by choosing covering hyperedges from the 
left to the right until all vertices in Conn are covered. The first separator returned 
by cover is {A,F}, the second one is {A, B}, etc. A complete procedure call 
sequence based on the above heuristics and the corresponding decomposition of 
the hypergraph in Fig. 1 is now shown in Fig. 6. 

In addition to this heuristic approach, we use also some kind of random- 
ization in our algorithm. In particular, the hyperedges read from the input file 
are randomly ordered. This guarantees that the results are independent from a 
particular representation of the hypergraph in the input file. However, this also 
means that det-fc-decomp has to be applied several times for evaluation pur- 
poses since the results and execution times may differ in each run. Of course, this 
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randomization step can also be removed to take the input order of hyperedges 
and vertices into account which sometimes improves the obtained results. 

5 Experimental Results 

In this section, we present our experimental results of the algorithm det-fc-de- 
comp when applied to three different classes of hypergraphs. We compare these 
results with the results obtained from evaluations of opt-fc-decomp (which can 
be downloaded from http://www.deis.unical.it/scarcello/Hypertrees/) 
and McMahan's heuristic approach. The latter one is based on Bucket Elimi- 
nation (BE) with ordering heuristics maximum- cardinality search, min-fill, and 
min-induced width [8]. Our evaluations were performed on two different ma- 
chines since opt-fc-decomp is only available as Microsoft Windows executable. 
So opt-fc-decomp was evaluated under Microsoft Windows 2000 on a 2.4GHz 
Intel Pentium 4 processor with 512MB main memory and both det-fc-decomp 
and BE were evaluated under SuSe Linux 9.2 on a 2.2GHz Intel Xeon processor 
(dual) with 2GB main memory. Note that the different memory sizes have no 
relevant influence on the results since the memory usage of det-fc-decomp and 
BE is far less than 512MB when applied to our examples. We chose the Intel 
Xeon as reference machine for normalizing the execution times. 

Since BE and det-fc-decomp use some kind of randomization, we applied 
them five times to each example in order to obtain representative results. From 
these five runs we selected the minimal width and computed the average run- 
time. Moreover, we defined a timeout of 3600 seconds for each run. Since we 
consider the width as the most important evaluation criterion and not the run- 
time, we set k always to the smallest value for which we obtained a solution from 
det-fc-decomp within our timeout. 

The first class of our benchmarks are from DaimlerChrysler and consist of 
hypergraphs extracted from adder and bridge circuits, the NewSystem examples, 
and a model of a jet propulsion system [6, 1] . The results when applying our three 
algorithms to theses hypergraphs are shown in Table 1. The first row in each line 
contains the name of the example and its size, i.e., the number of hyperedges 
(atoms) and the number of vertices (variables). The second row contains the 
hypertree- width if known, i.e., the optimal result. Then there are two rows for 
each algorithm; the first one contains the minimal width and the second one the 
average runtime in seconds. It is easy to see that det-fc-decomp outperforms 
both fc-decomp and BE on the DaimlerChrysler examples. 

The second class of benchmarks are hypergraphs extracted from two dimen- 
sional grids [1] . The advantage of these examples is that their hypertree- width is 
known by construction. Although we can easily construct an optimal decompo- 
sition of these examples by hand, it is seemingly very hard for our algorithms as 
can be seen in Table 2. Note that opt-fc-decomp cannot solve any of them within 
our timeout of one hour. Moreover, note that the runtime of det-fc-decomp is 
very high for larger examples. One reason for this may be that we set fc to 



15 



Instance (Atoms / Variables) 


Min 


opt-fc-decomp 


BE 


det-^-decomp 


Width 


Time 


Width 


Time 


Width 


Time 


adder .15 (76 / 106) 


2 


2 


2 


2 





2 





adder_25(126/ 176) 


2 


2 


20 


2 





2 





adder JO (251 /351) 


2 






2 





2 





adder_75(376/526) 


2 






2 





2 





adder 99 (496/694) 


2 






2 


1 


2 





bridgeJS (137 / 137) 


2 




9 


3 





2 





bridge_25 (227 / 227) 


2 


2 


69 


3 





2 





bridge_50 (452 / 452) 


2 


2 


1105 


3 


1 


2 





bridge_75 (677 / 677) 


2 






3 


1 


2 





bridge_99 (893 / 893) 


2 






3 


2 


2 


1 


NewSysteml (84/ 142) 


3 






3 





3 





NewSystem2 (200 / 345) 


3 






4 





3 





NewSystem3(278/474) 








5 


1 


4 





NewSystem4(418/718) 








5 


2 


4 





atv_partial_system (88 / 125) 


3 






3 





3 






Table 1. DaimlerChrysler benchmarks 



Instance (Atoms / Variables) 


Min 


opt-fc-decomp 


BE 


det-^-decomp 


Width 


Time 


Width 


Time 


Width 


Time 


grid2dJ0(50/50) 


4 






5 





4 





grid2dJ5(112/ 113) 


6 






8 





6 


3 


grid2d_20 (200 / 200) 


7 






12 





7 


3140 


grid2dJ5 (312/313) 


9 






15 


3 


10 


2000 


grid2d_30 (450 / 450) 


11 






19 


7 


13 


1566 


grid2dJ5 (612/613) 


12 






23 


15 


15 


1905 


grid2d_40 (800/800) 


14 






26 


28 


17 


2530 


grid2d_45 (1012/ 1013) 


16 






31 


51 


21 


2606 


grid2d_50 (1250/ 1250) 


17 






33 


86 


24 


2786 


grid2d_60 (1800/ 1800) 


21 






41 


204 


31 


2984 


grid2d_70 (2450/2450) 


24 






48 


474 


42 


2161 


grid2d_75 (2812/2813) 


26 






48 


631 


45 


2881 



Table 2. Grid2D benchmarks 
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the smallest value for which we obtain a decomposition within an hour. If we 
increase k, then the runtime will decrease in many cases. 



Instance (Atoms / Variables) 


Min 


opt-k-decomp 


BE 


det-A>decomp 


Width 


Time 


Width 




Width 


Time 


cT7 (\ 1 1 17") 
SZ / (13/ ll) 
















cTOS ( 1 C\A 1 1 1 ^\ 
SZUo(lU'+/ lij) 


— 






7 





6 





,.1110 (" 1 / 1 t(i\ 


> 3 














WW { I / J / 1 OH- _) 


— 














c ( 1 7fi / 1 S^l 


— 














cIRT ( 1 7Q / 1 ST\ 

s^oz (i /y / lozj 


— 














sJ>oO ( 1 uj / 1 / ZJ 








8 


1 


7 


1824 


s400(183/ 186) 


> 3 






6 





5 


273 


s420 (212/231) 


> 3 






9 





8 


454 


s444(202/205) 


> 3 






6 





5 


385 


s510(217/236) 


> 3 






23 


1 


20 


2082 


s526 (214/217) 


> 3 






8 


1 


7 


1715 


s641 (398 / 433) 








7 


1 


7 


1611 


s713(412/447) 








7 


1 


7 


1800 


s820 (294/312) 


> 3 






13 


3 


12 


2846 


s832(292/310) 


> 3 






12 


3 


11 


2575 


s838(422/457) 


>3 






16 


1 


15 


2046 


s953(424 / 440) 


> 3 






40 


8 






S1196 (547/561) 








35 


11 






S1238 (526 / 540) 








34 


13 






S1423 (731/748) 








18 


3 






S1488 (659/667) 








23 


18 






sl494(653/661) 








24 


19 






S5378 (2958 / 2993) 








85 


141 







Table 3. ISCAS89 benchmarks 



Finally, the third class of benchmarks arc hypcrgraphs extracted from circuits 
of the ISCAS89 (International Symposium on Circuits and Systems) benchmark 
suite. The ISCAS benchmarks are again examples from practice. However, they 
are obviously much more difficult to solve than the DaimlcrChrysler examples 
as can be seen in Table 3. Note that the widths obtained from det-fc-decomp 
are not much smaller than those from BE although it needs much more time. 

In summary, we can conclude that det-fc-decomp significantly outperforms 
opt- fc-decomp on all classes of benchmarks we have tested, opt-fc-decomp ter- 
minates within our timeout only for very small and simple examples. In contrast, 
det-fc-decomp terminates also for larger examples with results comparable to or 
even better than those computed by BE. Thus, if the width is the most important 
criterion and not the runtime, det-fc-decomp also outperforms BE as long as the 
hypergraphs are not too large and too complicated. 
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6 Conclusion 



We have presented the new algorithm det-fc-decomp for constructing hyper- 
tree decompositions of hypcrgraphs. This algorithm results from the alternating 
algorithm fc-decomp introduced by Gottlob et al. [3] when replacing the guess 
and c/iecfc-part by a search procedure. We have evaluated det-fc-decomp by a 
series of benchmark examples from industry and academics and compared it to 
the currently best hypertree decomposition algorithms. Our results have shown 
that det-fc-decomp performs very well; in particular, it significantly outper- 
forms the algorithm opt-fc-decomp. Future work is to further restrict the search 
space by theoretical results and to improve the heuristics used to accelerate 
the search process. 
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A APPENDIX 



Proposition 1. The algorithm fc-decomp presented in Algorithms 1 and 2 is a 
notational variant of the algorithm fc-decomp presented in [3]. 

Proof. First note that fc-decomp returns NULL if it rejects and it returns a c- width 
hypertree decomposition if it accepts, where c < fc. Algorithm 1 represents the 
main procedure which calls the recursive procedure fc-decomposable described 
in Algorithm 2. The names of these procedures are the same as in their original 
presentation in [3]. The first parameter of fc-decomposable, however, is a set 
of edges instead of a set of vertices. In particular, we use the set edges (Cr) of 
hyperedges containing some vertices in the component Cr instead of the compo- 
nent Cr itself. Both variants are equivalent since Cr = (J edges (Cr ) \ [J R can 
be computed from edges{CR) and R. So the first parameter of fc-decomposable 
is the set of hyperedges representing the current component and the second pa- 
rameter is the set of A-labels of the parent node. Note that we use the term 
separator for a set of hyperedges chosen as A-labels since they separate the sub- 
components. 

Now, let us compare Algorithm 2 with its original formulation in [3]. The 
guess-step of Separator in line 1 is obviously the same. For the check-step in 
lines 2 to 4, it is easy to see that both conditions are equivalent to their orig- 
inal formulation. Line 6 describes the computation of the set C of components 
in [3]. This computation is now hidden in the procedure separate and its re- 
sult Components is a set of sets of hyperedges such that for each C E C there 
exists Comp E Components such that \J Comp \ (J Separator = C. The remainder 
of Algorithm 2 contains the recursive call as in [3] and the construction of the 
hypertree decomposition. In particular, the recursive call itself is given in line 9, 
where we loop through all components. If the decomposition of a component is 
rejected, we reject in line 11; otherwise, we store the hypertree decomposition 
of the component in Subtrees in line 13. Finally, a new hypertree node is con- 
structed in line 17 with Separator as A-labels and Subtrees as subtrees; the 
X~ labels are computed according to line 16. Note that the \~ labels computed in 
this way do not necessarily satisfy the third normal-form condition. However, 
this condition can be easily satisfied by a simple post-processing step. 

So we have demonstrated the equivalence between our notational variant of 
the algorithm fc-decomp and its original formulation in [3]. □ 
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